Burst statistics in Alcator C-Mod SOL turbulence 
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Bursty fluctuations in the scrape-off layer (SOL) of Alcator C-Mod have been analyzed 
£T) \ using gas puff imaging data. This reveals many of the same fluctuation properties as Lang- 

muir probe measurements, including normal distributed fluctuations in the near SOL region 
while the far SOL plasma is dominated by large amplitude bursts due to radial motion of 



blob-like structures. Conditional averaging reveals burst wave forms with a fast rise and 
slow decay and exponentially distributed waiting times. Based on this, a stochastic model 
of burst dynamics is constructed. The model predicts that fluctuation amplitudes should fol- 

'i-H . 

low a Gamma distribution. This is shown to be a good description of the gas puff imaging 
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data, validating this aspect of the model. 
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Q\ I. INTRODUCTION 

\o 

iT) \ Cross-field transport of particles and heat in the scrape-off layer (SOL) of magnetically con- 

O 

fined plasmas is dominated by radial motion of blob-like structures QJJ— 13Q - The average cross-field 
particle and heat fluxes caused by such filaments depend on their amplitude distribution and fre- 
^ quency of occurrence. The statistical properties of plasma fluctuations in the SOL are thus crucial 

for development of a first-principles physics based description of transport and main-chamber in- 
teractions, and may also prove important for understanding the empirical discharge density limit 

Im. 

Plasma fluctuations in the Alcator C-Mod SOL have been investigated by analysis of data from 
gas puff imaging (GPI) measurements at the outboard mid-plane region in a set of ohmically 
heated, lower single null discharges with a scan in line-averaged density. It is shown that the 
GPI diagnostic reveals many of the same fluctuation properties as Langmuir probe measurements, 
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including normal distributed fluctuations in the near SOL region while the far SOL plasma with 
a broad particle density profile is dominated by large amplitude bursts due to radially outward 
motion of blobs IjJ-y]. The fluctuation probability distributions in the far SOL are strongly skewed 
and flattened with an exponential tail towards large values for all line-averaged particle densities 

uUao. 

The burst statistics are revealed by means of a standard conditional averaging technique. GPI 
measurements reproduce the familiar asymmetric burst wave forms with a fast rise and a slow de- 
cay and an average duration that is independent of the line-averaged density [|2j]. The waiting times 
between large amplitude events are found to be exponentially distributed for all line-averaged den- 
sities. Based on these results, a stochastic model for the intermittent fluctuations has been con- 
structed. This model reveals the importance of burst duration, waiting times and amplitudes for 
large far SOL plasma densities and fluctuation levels. From the model it is shown that the fluctua- 
tion amplitudes follow a Gamma distribution and that there is a parabolic relationship between the 
skewness and kurtosis moments. This compares favourably with the GPI measurements. 



H. EXPERIMENTAL SETUP 

For this investigation we analyze GPI data from a set of similar, deuterium fuelled, ohmically 
heated, lower single null plasmas with plasma current 800 kA and toroidal magnetic field 4.0T. 
The GPI diagnostic consists of a 9 x 10 array of in- vessel optical fibres with toroidally viewing, 
horizontal lines of sight which are locally enhanced in the object plane by an extended He gas puff 
from a nearby nozzle. Excitation of the He neutral gas, and thus the intensity of the GPI signals, 
is determined by a combination of the local electron particle density and temperature [4]. The 
fibres are coupled to high sensitivity avalanche photo diodes and the signals are digitized at a rate 
of 2 x 10 6 frames per second. The viewing area covers the major radius from 88.00 to 91.08cm 
and vertical coordinate from —4.51 to —1.08 cm with an in-focus spot size of 3.7 mm for each of 
the 90 individual channels. The radial position of the last closed flux surface at the vertical centre 
of the image, Z = —2.61 cm, is in the range from 89.4 to 89.7 cm for all the discharges presented 
here. The limiter radius mapped to this vertical position is at R = 91.0cm. Further information 
about the GPI diagnostic can be found in Ref. |4|]. 

This paper reports on experiments performed during the FY20 10 run campaign using GPI data 
from run 1 100803. This comprises a four-point scan in line-averaged density with the Greenwald 
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fraction n e /«G from 0. 15 to 0.30. Here the Greenwald density is given by hq — I v / (na 2 ) 10 20 m~ 3 
where I p is the plasma current in units of MA and a is the plasma minor radius in units of meters. 
For the present density scan hq = 5.26 x 10 20 m~ 3 . The condition at the outer divertor goes from 
sheath limited at the lowest density to high recycling at the highest density in this scan. For 
each discharge the GPI diagnostic yields 0.25 s usable data time series during the flat-top of the 
plasma current. By combining data from two discharges at the same n e and two nearby diode 
channels with identical statistical properties, we obtain single-point time series of one second 
duration which allows calculation of statistical averages with high accuracy. 



HI. GPI MEASUREMENTS 

In Fig. Q] we show the radial variation of the relative fluctuation level of the GPI intensity signals 
at Z = —2.61 cm. The fluctuations increase drastically in magnitude with radial distance into the 
SOL. As shown in Fig. [2l the raw time series are here dominated by large- amplitude bursts due 
to the radial motion of blob-like structures Q5J]. This results in positively skewed and flattened 
probability density functions (PDFs) of the intensity signals. As an example of this, we present in 
Fig.[3]the distribution function for the GPI signals for n e /nQ = 0.20 at Z = —2.61 cm for the four 
GPI fibre view positions radially outside the separatrix. It is clearly seen that the PDF changes 
from a normal distribution in the near SOL region to strongly skewed and flattened in the far SOL. 
At the limiter radius R = 91 cm, the PDF has an exponential tail towards large signal amplitudes. 
This is similar to what has previously been found from Langmuir probe measurements []2l |3||. 

The radial motion of blob-like structures through the SOL results in single-point recordings 
dominated by bursts with a fast rise and a slow decay. This is demonstrated by the asymmetric 
wave form obtained from conditional averaging presented in Fig.[4j Here it is seen that the average 
burst duration is the same for all line-averaged densities. This is again similar to what has pre- 
viously been found from probe measurements as well as numerical turbulence simulations D2j,|6J]. 
The waiting time between large amplitude bursts is also obtained from the conditional averaging. 
As shown in Fig. |51 the burst waiting times are found to be exponentially distributed for all line- 
averaged densities. An exponential distribution describes the time between events for a Poisson 
process, in which events occur randomly and at a constant average rate. The large-amplitude blobs 
appearing in the far SOL are thus uncorrelated. 
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IV. STOCHASTIC MODELLING 

The experimental measurements presented here suggest that fluctuations in the SOL can be 
represented as a random sequence of bursts similar to the classical "shot noise" process D7|], 

®(t)=Y i AkY(t-tk), (1) 

k 

where tk is the burst arrival time for event k and i// is a fixed burst wave form. The burst waiting 
times, given by Tk = h~ h-u are assumed to be exponentially distributed, consistent with the 
results presented in Fig. |5] The individual burst wave forms will for simplicity be approximated 
by a sharp rise followed by a slow exponential decay, y(t) = 0(f) exp(— t/x^), where is the step 
function and the burst duration time Td is taken to be constant. The burst amplitudes A are also 
taken to be exponentially distributed, 

P A (A) = -i- exp ( ~\ (2) 



(A) r V (A) 

where (A) is the average burst amplitude. For the exponential burst wave form, the mean value of 
the signal in Eq. (Q} is readily calculated 0,180, 

(<£> = ^<A>, (3) 

where T w is the average burst waiting time. The above equation elucidates the role of burst dura- 
tion, waiting time and amplitude for large SOL plasma densities. 

In the case of exponentially distributed burst amplitudes, the relative fluctuation level can be 
written as ^nns/ = ( ivy/id) 1 / 2 , while the skewness and kurtosis moments for 4> are given by 

u 



4t w V /2 _ - 6t, 



5=1 — 1 , K = 3 + ^. (4) 

This shows that there is a parabolic relation between these moments given by K = 3 + 3S 2 /2. Like 
the relative fluctuation level, also the skewness and kurtosis increases with the ratio T w / Td- The 
parameter y = x&/ T w is thus a measure of intermittency in the signal given by Eq. ©. 

In Fig. [6] we present the kurtosis as function of skewness calculated for all GPI fibre view 
positions located in the SOL for the discharges in the present density scan. It is clearly seen that 
the skewness and kurtosis moments increase radially outwards in the SOL, similar to the relative 
fluctuation level presented in Fig. \T\ Also shown in Fig. |6]is the predicted parabolic relation given 
by Eq. ©. This is an excellent description of the experimental data for all but the outermost 
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diode channel positions. Reduced emission from the neutral He gas is expected for the cold and 
dilute plasma in the limiter shadow. Here significant emission levels only arise in the presence 
of large amplitude blob structures, which leads to excessively large higher order moments of the 
fluctuations. This is likely the cause of the apparent discrepancy from a parabolic relation for the 
signals from the outermost GPI views in Fig. |6] Such a parabolic relation between the skewness 
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and kurtosis moments has been found for many other plasma experiments 

For exponentially distributed burst waiting times and amplitudes, the PDF for 4> given by the 
model (OQ) can be shown to be a Gamma distribution Q8D , 

where the scale parameter is given by (<I>) /y and the shape parameter is y = (4>) 2 /4>^ ris . In Fig. [3] 
we have presented the corresponding Gamma distribution, which is seen to be an excellent fit to 
the data for the innermost SOL positions. Note that the PDFs for the experimental data comprise 
more than four decades on the ordinate axis — a consequence of the long data time series used in 
this analysis. In the limiter shadow the PDF of the measured signals have a strongly elevated tail 
which is likely due to suppressed emission in the cold far SOL plasma as discussed above. The 
Gamma distribution has previously been found to accurately describe the ion saturation current 
signal in the SOL of TCV across a broad range of plasma parameters [|1( 



V. CONCLUSIONS 



GPI measurements in the SOL of Alcator C-Mod have shown that plasma fluctuations are 
dominated by large amplitude bursts due to radial motion of blob-like structures. The burst wave 
form is asymmetric with a fast rise and slow decay. Conditional averaging reveals that the burst 
waiting times and amplitudes are exponentially distributed, thus blobs appearing in the far SOL 
are uncorrelated. The PDF of the GPI intensity signal changes from a normal distribution in the 
near SOL to strongly skewed and flattened with an exponential tail in the far SOL. A stochastic 
model of burst dynamics is constructed based on three parameters: burst duration, waiting time and 
amplitude. Consistent with the predictions of the model, the PDFs of the GPI intensity fluctuations 
are well described by a Gamma distribution and there is correspondingly a parabolic relation 
between the skewness and flatness moments. These results indicate that the stochastic model is an 
accurate description of the burst dynamics in the tokamak SOL. 
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FIG. 1. Radial variation of relative GPI intensity fluctuation level at Z = —2.61cm. The position of the last 
closed magnetic flux surface is indicated by the shaded region. 




FIG. 2. Intensity signal for n e /n G = 0.20 at (R,Z) = (90.69, -2.61) cm. 




FIG. 3. PDF of GPI intensity signals for n e /nc = 0.20 at Z = —2.61cm and various radial positions in the 
SOL. The solid lines are predictions of the model. 



8 




FIG. 4. Conditionally averaged wave form for GPI intensity burst amplitudes larger than 2.5 times the mean 
value at the reference position (R,Z) = (90.69, —2.61) cm. 
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FIG. 5. PDF of waiting times between burst events with maximum amplitude larger than 2.5 times the mean 
value at the reference position (R,Z) = (90.69, —2.61) cm. 
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FIG. 6. Kurtosis versus skewness for the GPI intensity signals in the SOL for all discharges in the density 
scan. The full line is the prediction of the model. 



